Kerr black hole lensing for generic observers in the strong deflection limit 
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We generalize our previous work on gravitational lensing by a Kerr black hole in the strong 
deflection limit, removing the restriction to observers on the equatorial plane. Starting from the 
Schwarzschild solution and adding corrections up to the second order in the black hole spin, we 
perform a complete analytical study of the lens equation for relativistic images created by photons 
passing very close to a Kerr black hole. We find out that, to the lowest order, all observables 
(including shape and shift of the black hole shadow, caustic drift and size, images position and 
magnification) depend on the projection of the spin on a plane orthogonal to the line of sight. In 
order to break the degeneracy between the black hole spin and its inclination relative to the observer, 
it is necessary to push the expansion to higher orders. In terms of future VLBI observations, this 
implies that very accurate measures are needed to determine these two parameters separately. 

PACS numbers: 95.30.Sf, 04.70.Bw, 98.62.Sb 
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I. INTRODUCTION 



As predicted by General Relativity, photons passing 
near a black hole suffer deviations from their original 
trajectory. If the minimum distance between photon and 
black hole is much larger than the gravitational radius, 
a weak field approximation of the metric tensor is suffi- 
cient to describe the light deflection. Two images of the 
original source are then detected by the observer. On 
the other hand, photons passing very close to the black 
hole may suffer very large deviations without falling into 
the black hole. These photons may perform one or more 
loops around the black hole before reemerging in the ob- 
server direction, thus generating two infinite sets of rel- 
ativistic images very close to the black hole shadow. It 
can be easily intuited that these relativistic images rep- 
resent a unique probe to gain information on the very 
strong gravitational fields surrounding the black holes. 
Through their study it would be possible to learn the 
properties of black holes and get new insight on General 
Relativity in a strong field regime. The features of rel- 
ativistic images will thus represent a possible challenge 
arena for alternative theories of gravitation. 

Even though a general relativity approach to this sub- 
ject typically results in involved equations and heavy nu- 
merical integrations, a surprisingly simple formula for the 
deflection angle induced by a Schwarzschild black hole 
in the Strong Deflection Limit (SDL) was proposed by 
Darwin jj] and revived in Refs. [1, 0, |j|- The loga- 
rithmic divergence of the deflection angle in the impact 
parameter was showed to be universal and valid for all 
spherically symmetric black holes Q, as every class of 
such black holes leads to the same expansion for the de- 
flection angle, with coefhcients depending on the specific 
form of the black hole metric. The SDL method was 
then applied to several classes of black holes, ranging 



from Reissner-Nordstrom to black holes in string theory, 
from braneworld black holes to wormholes Q. By the 
SDL method it is thus possible to quantify the observ- 
ables related to relativistic images for any class of spher- 
ically symmetric black holes, allowing an easy compar- 
ison among different theories. For alternative methods, 
see Refs. 7, 8]. 

For spinning black holes, things do not work so eas- 
ily. Starting from the geodesies equations in Kerr space- 
time, that Carter [2j reduced to first order equations de- 
pending on four constants of motion, many numerical 
approaches have been developed to study and visualize 
such geodesies. Numerical efforts have also been profused 
in the context of gravitational lensing to investigate the 
apparent shape of the accretion disk of the black hole 
[ 4 lift llfj . the light curve of a star orbiting around it 
[l2 | and the structure of the caustics [Ij| , which turned 
out to be extended and to have a 4-cusped astroid struc- 
ture. Some interesting general results have recently been 
derived through Morse theory [lj| . The extension of the 
SDL methodology to Kerr black holes was firstly per- 
formed in Ref. 15] and the SDL formula was recovered 
for photons lying near to the equatorial plane. Anyway 
the expansion coefficients had to be calculated numeri- 
cally as functions of the lens spin. 

A first step toward a complete analytical treatment of 
this subject was made in Ref. [13 (hereafter Paper I) 
where the lens equation was analytically solved in the 
limit of small values of the angular momentum of the 
black hole (denoted by a) and for observers lying on its 
equatorial plane. This last assumption, besides ensuring 
simpler equations, was justified by the fact that the most 
important black hole candidate (Sgr A*, firstly suggested 
in Ref. iq|) is hosted in the center of our Galaxy and 
presumably has a spin-axis perpendicular to the Galactic 
plane, where the solar system lies. The expansions for 
small values of the angular momentum allowed to use 



the Schwarzschild SDL formula as a starting point for the 
description of the deflection of hght rays looping around 
a Kerr black hole. 

This analytical approach provided very simple equa- 
tions (which could even be inverted for sources near to 
a caustic) and a full description of the extended struc- 
ture of the caustics, which were confirmed to have a 4- 
cusped structure, symmetric w.r.t. the equatorial plane 
and shifted from the optical axis. Only the first order 
caustic cannot be recovered in the SDL approach as it 
is formed in the weak deflection regime [liil \ll\ . It was 
also shown that the extension of relativistic caustics en- 
hances the cross section for the creation of additional im- 
ages, whose magnification is sensible in a relatively large 
region around the caustic. Direct observations of these 
relativistic images, which should be possible with the res- 
olutions achieved by future projects, could test the Kerr 
nature of black hole lenses (see e.g. Refs. [l^ [l^ for 
detailed discussions on observational perspectives). It is 
interesting to compare the lensing effect of a Kerr black 
hole to that of a Schwarzschild black hole embedded in an 
external gravitational field. Also in the latter case astroid 
caustics arise, though with different sizes and positions 

In this Paper we further investigate Kerr black hole 
lensing, getting rid of the equatorial ofeseruer hypothesis. 
In spite of the presence of a new parameter (the incli- 
nation of the spin axis relative to the line of sight, that 
we shall indicate by "(Jo), the surprisingly simple struc- 
ture of all analytical results is preserved. Our philosophy 
will be to try to confine all technicalities to the appen- 
dices or refer the reader to Paper I for more detailed 
derivations. This paper will thus keep its main focus on 
the implications of all results on observable quantities. 
What emerges from our study is that all observables (to 
the lowest order) just depend on asin^o, that is the pro- 
jection of the spin on a plane orthogonal to the line of 
sight. The consequences of this fact will be discussed in 
the conclusions in Section 7. 

Our paper is structured as follows: in Section 2 we 
recall the main properties of Kerr geodesies. In Section 
3, we trace the borders of the shadow of the Kerr black 
hole for all values of the observer declination. In Section 
4 we apply the SDL to null Kerr geodesies illustrating 
the main strategy and referring to two appendices for 
the details. In Section 5 we derive the critical curves and 
caustics structure and in Section 6 we analyze the lens 
mapping in the neighbourhood of a caustic, finding the 
position and the magnification of the images, conclud- 
ing with a discussion on the detectability of relativistic 
images. 



II. KERR GEODESICS 

In this section, we shall review the basics of Kerr 
geodesies and introduce the notations to be used through- 
out the paper. For more details on the physical meaning 




FIG. 1: Boyer-Lindquist coordinates in a Kerr metric, also 
referred as spin-oriented coordinates in the text. L is the 
black hole with spin a. O is the observer and P is a generic 
point. The gray disk visualizes the equatorial plane of the 
black hole. 



of all quantities, the reader may refer to Paper I or di- 
rectly to Ref. [23. 

The main subject of our paper is the Kerr black hole, 
whose metric in Boyer-Lindquist coordinates J2j|, a;'^ = 
(t, X, d, 4>) reads 



ds^ = dt^ - '--dx^ - p^dd^ 



p'- A 

(:E2 + a2)'-a2Asin2,? 



2 „a jj,2 



sin ■dd(l) 



lax sin t? 



dtd(f> 



A= x^-x + a"^, 
p^ = x^ + a^ cos^ -d. 



(1) 

(2) 
(3) 



Distances are measured in units of the Schwarzschild ra- 
dius {2MG/c^ = 1), 1? and <p are the colatitude and az- 
imuth respectively, x is the radial coordinate and a is 
the specific angular momentum of the black hole, run- 
ning from (Schwarzschild black hole) to 1/2 (extremal 
Kerr black hole) in our units. 

We consider a static observer at Boyer-Lindquist co- 
ordinates {DoL,'&o,4>o)- The distance between observer 
and black hole is thus Dql, while the colatitude ^o of 
the observer coincides with the inclination of the spin on 
the line of sight OL. Exploiting the freedom to choose 
the zero of the azimuth, we set (f)o = tt. We will very 
often use the notation /i = cosi5. Thus we also define 
Po = cos -do- Fig. n illustrates Boyer-Lindquist coordi- 
nates for a generic point P and for the observer O in 
particular. 

Lightlike geodesies can be expressed in the following 
form in terms of the first integrals of motion J and Q 



found by Carter Q 
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and 4>i is the initial value of the azimuthal coordinate of 
the photon. 

The roots of R represent inversion points in the ra- 
dial motion. In gravitational lensing we consider photons 
coming from infinity, grazing the black hole and going 
back to infinity. For such trajectories there is only one 
inversion point a;o , representing the closest approach dis- 
tance. The minimum allowed value of xq can be found 
solving the equations R{x) = and R'{x) — simulta- 
neously. However, in Kerr black hole, we do not have a 
unique minimum closest approach Xm , but rather a con- 
tinuous family of values which depend on the approach 
trajectory followed by the photon. In particular, it is pos- 
sible to establish a relation among the minimum closest 
approach Xm and the corresponding values of the con- 
stants of motion J and Q, that we shall indicate by Jm 
and Qm (see e.g. Ref. 20|) 






X^,^{2Xm-3)+a^{l + 2Xrn) 

a{l - 2xm) 

xl^[2a^~Xrn{Xm-S/2)^] 
a^(Xra-\l2f 



(8) 
(9) 



Xm also represents the radius of the unstable circular 
photon orbit. This radius is fixed to 3/2 when a = 
(Schwarzschild black hole). In the case of Kerr black 
holes, Xra may vary between two limiting values Xm+-, 
Xm-, depending on the incoming direction of the photon. 
The two limiting values can be analytically obtained solv- 
ing the equation Qm = (in fact, it is possible to show 
that gravitational lensing trajectories cannot be realized 
for Q < ^). To the third order in a, they read ^ 
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+ 0(a*). (10) 



For example, photons whose orbit lies on the equato- 
rial plane may turn either in the same sense of the black 
hole (prograde photons) or in the opposite sense (ret- 
rograde photons). Prograde photons are allowed to get 
closer to the black hole, with a minimum closest approach 
given by Xm+, while retrograde photons must stay far- 
ther than Xm-, in order to be defiected without falling 




FIG. 2: The shadow of the black hole in the observer sky for 
a — 0.1 and different values of the observer position i?o. The 
solid line is for i9o ~ 7r/2 (equatorial observer), the dashed 
line is for i9o = tt/I and the dotted line is for i?o = (polar 
observer) . 



into the black hole. Photons whose orbit does not lie 
on the equatorial plane are characterized by intermedi- 
ate values of Xm, with Qm > 0. Thus Xm can be used to 
parametrize the family of unstable photon orbits allowed 
in Kerr metric or, equivalently, the incoming direction of 
the photon. The corresponding values of the constants 
of motion are uniquely determined by Eqs. ||SJ| and ©. 

Although exact expressions for Xm+ and Xm- are avail- 
able, it is convenient to start with a perturbative expan- 
sion ab initio in order to be prepared to face more com- 
plicated quantities in the following lo| . Throughout our 
treatment, only for Xm we need to push the expansion to 
the third order, in order to obtain some quantities to the 
second order in a. 



III. THE SHADOW OF A KERR BLACK HOLE 

The constants of motion J and Q have an immediate 
link to the position in the sky where the observer detects 
the photon. In fact, we can define angular coordinates 
(6*1, 6*2) on the observer sky centered on the black hole 
position. We choose the orientation of these coordinates 
in such a way that the spin axis of the black hole is pro- 
jected on the 6'2-axis (see Fig|2)l. 

As shown in Ref. [23], photons detected by the ob- 
server at angular coordinates (^1,^2) are characterized 



by constants of motion given by 



J 



-eiDoLV^-fj-l 



Q = elDl^ + ^Jil{elDl^~a'). 



(11) 

(12) 



These relations can be easily recovered taking the limit 
for large distances in the equations of motion of the pho- 
ton. They show that J can be identified with the com- 
ponent of the orbital angular momentum of the photon 
along the spin axis, whereas Q + J'^ + ii^a? is the squared 
total angular momentum of the photon. 

Note that, with our choice of (6*1,02), in the limit of 
equatorial observer /Xo = 0, prograde photons (J > 0, 
(5 = 0) are detected by the observer on the left side of 
the black hole, while retrograde photons (J < 0, Q = 0) 
are detected on the right side. Conversely, in the limit 
of polar observers {^.q — + ±1), the projected angular mo- 
mentum J vanishes, while Q ^r {Ol -\- 02) ^ol ~ '^ ■ 

Inverting Eqs. pi(l and H12|l . wc find the position 
(^1, 6*2) in the sky where the photon is detected with given 
constants of motion J and Q, apart from an ambiguity 
in the sign of Qi 



J 



DolV^^I^o' 
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(13) 
(14) 



These relations can be used to convert the locus 
[JrmQm), parametrized by Xm according to Eqs. © 
and jni in the ( J, (3)-space, into a new one (0i,m,02,m) 
in the observer sky. However, not all values of Xm in the 
range [xm+TXm-\ are acceptable. This can be easily un- 
derstood, as photons lying on the equatorial plane can 
never reach non-equatorial observers. The reality condi- 
tion for 02, m restricts Xm to the range [a;p+,a;j,_], where 
Xp+ and Xp- are the roots of the equation 6*2,™ = 0. To 
third order in a, these quantities read 



Xp±= - T -Jf^^ - A^o - ga^(l + Mo) 



T 



27^3 



\b + Q^il)^T^, + 0{a^). (15) 



Comparing with Eq. pOfl . we see that Xp± — > Xm± in 
the limit /io — >■ 0. On the other hand, when /Xq — > ±1, 
the allowed range for Xm shrinks to a single value Xp -^ 
I — |a^. This witnesses that when the observer is on the 
polar axis the axial symmetry of the lensing configuration 
is restored and all unstable photon orbits have the same 
radius again. 

When a vanishes, Xp+ and Xp_ both coincide with the 
Schwarzschild photon sphere radius, 3/2, while, when a 
is not zero, they are distinct and every value of Xm in the 
interval [xp+,a;p_] uniquely fixes the amplitude of the 
oscillation of the photon orbit on the equatorial plane, as 
we shall see later. On the basis of this consideration, in 
Paper I (with ^lq = 0) we introduced a parametrization 



of Xra in the range [Xm+, Xm-\^ replacing a with aS, in Eq. 
(|10|l . with the parameter ^ varying in the range [—1, 1]. 
In order to take into account the changes from Eq. H10|l 
to (|15|) . we have to update such parametrization, since it 
is not directly applicable to the case /Xq 7^ 0. Our new 
parametrization for Xm is 



2-^<yT^-9«^(i 



^^l) 



27\/3 



'a^ + Q^il)^T^, + 0(a^). (16) 



As ^ varies in the interval [—1,1] we get all possible 
values of Xm in the interval [xp+,XpJ\. It will become 
clear later that ^ is strictly related to the position angle 
of the generic point in the observer sky. 

With this parametrization, we can rewrite Eqs.©-® 
to the second order in a as 
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Qm(0 
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(18) 



Notice that the presence of a in the denominators of 
Eqs. 10)- (O allows ^ to be present already in the zero- 
order terms in Eas. l|17|) - (|18|l . permitting the use of the 
^-parametrization in Schwarzschild spacetime as well. 
However, since this parametrization has been introduced 
in a slightly different way w.r.t. Paper I, the expressions 
derived here cannot be directly compared to those of Pa- 
per I, except for those quantities that are independent 
of ^. For example, eliminating ^ from Eqs. (|17|l and 
(fT5|l . one can derive an expression for the locus (J„j, Qm) 
in the form Qm{Jm)- Doing the same with the expres- 
sions of Paper I, one would indeed find the same function 
Qm{Jin) in the limit ^o -^ 0. 

Inserting Eqs. ^\ and I^HJ in Eqs. I^SJ-ClIl we get 
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[l + 2^ll-2e{l-^il)]. (20) 



This locus is formed by the points in the observer sky 
where photons with minimum closest approach would 



show up. No gravitational lensing images are possible 
inside this locus, which is thus also known as the shadow 
of the black hole. In Fig. El we show it for different 
values of ^o- Note that, to zero order, 9i^m esc — ^ and 
02,m oc \/l — ^^, justifying the identification of ^ with the 
cosine of the position angle in the (6'i, 6*2) plane as taken 
from the opposite of the 6'1-axis. This fact facilitates the 
physical interpretation of the parameter ^. 

The shadow of the black hole is the first observable 
in extreme gravitational lensing by supermassive black 
holes. It thus deserves some further analysis in order to 
understand the effect of the spin and the observer posi- 
tion. 

First we note that 9i,m and 6*2,771, to second order in a, 
satisiy the ellipse equation 
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with the origin shifted rightward by 
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(21) 



(22) 



and semiaxes given by 



hole to such accuracy that we are able to extract a from 
a measure of the minor semi-axis Ai solely. However, 
since the spin contribution to the major semi-axis is only 
of second order in a, we need a very high accuracy in the 
shadow observation in order to appreciate such a small 
contribution. For example, if a = 0.1, the spin contribu- 
tion to Ai is of order 0.2%. As already pointed in Ref. 
|2J| by numerical examples, the disentanglement of a and 
■do is only possible for values of the black hole spin very 
close to the extremal case. By our perturbative formulae, 
we have justified this claim analytically. Of course, for 
high values of a, when higher orders contribute to deter- 
mine the shape of the shadow, the degeneracy between a 
and do can be broken, in agreement with what stated in 
Ref. "23. 

It has been pointed out in Paper I that as long as we 
deal with Kerr black holes with spin smaller than a — 
0.2, the perturbative approximation works surprisingly 
well. Then, the degeneracy between a and do in the 
shadow of the black hole poses a serious difficulty to the 
determination of the parameters of the black hole by the 
simple observation of the shadow. As we shall see in the 
next sections, this degeneracy plagues all gravitational 
lensing observables in different degrees. 
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By these analytical expressions for the shadow, we can 
make several considerations. The presence of a non- 
vanishing spin causes a slight distortion and a displace- 
ment of the shadow from the black hole position. When 
the observer lies on the spin axis {fio = il)i the axial 
symmetry is restored and the shadow returns to be cen- 
tered on the black hole and circular. However, even in 
this limiting case, the radius of the shadow is no longer 
3\/3/2 as in Schwarzschild, but it is slightly smaller, be- 
ing 3V3/2 - aVVS. 

It has been proposed that the observation of the shape 
of the shadow of a black hole by VLBI may help to de- 
termine the parameters of a Kerr black hole, such as its 
mass , its angular momentum and the inclination of the 
spin [lOl 12^. However, both in the shift Oq and in the 
ellipticity 
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the black hole spin and the observer declination appear in 



the same combination a^/l 



M^ 



a sin do, which repre- 



sents the projection of the spin on a plane orthogonal to 
the line of sight. Thus it is impossible to determine both 
the absolute value of the spin and its inclination from the 
shape of the shadow. The only possibility is that we al- 
ready know the distance Dql and the mass of the black 



IV. KERR LENSING IN THE STRONG 
DEFLECTION LIMIT 

As in Paper I, we introduce the following parametriza- 
tion of the observer sky 



92(e,0 



9i,m(C)(l 
92,™ (0(1 






(26) 



Varying ^ in the range [—1,1] and e in the range [—1, 00], 
we can obviously cover the whole upper half of the ob- 
server sky, since £, establishes the position angle of the 
light ray w.r.t. the (— 0i)— axis (through Eqs. l(T^ 
and H20I) ') whereas e fixes the angular distance from the 
shadow of the black hole. In fact, in the following, e 
will be generically referred to as the separation from the 
shadow of the black hole. 

We are interested into light rays experiencing very 
large defiections by a Kerr black hole. These light rays 
reach the observer from directions (6*1, 02) very close to 
the shadow. In the parametrization (|26|l . they are thus 
characterized by a very small positive e, while keeping 
^ in the whole range [—1, 1]. The SDL amounts to per- 
forming the integrals in the geodesies equations Q-lO, 
to the lowest orders in the separation from the shadow e. 

The values of the constants of motion J and Q, cor- 
responding to such strongly deflected photons, can be 
found using Eqs. l(TT|) - ((n|) : 






Q™(e)(l + 2e) 



2ahfil. 



(27) 
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Substituting these expressions in Eq. Q and solving 
the equation R = for xq, we get the closest approach 



distance as 

2 ; 2 



-/l 
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(29) 
(30) 



where we have introduced the compact notation 

i = Cv/t^m!. (31) 



As e represents the separation of the image in the ob- 
server sky from the shadow of the black hole, S represents 
the relative difference between the closest approach xq 
and the minimum closest approach Xm(0 fixed by the 
position angle through ^. It will be synthetically called 
approach parameter. As S decreases, we expect the de- 
flection to increase more and more. In the limit S —^ 0, 
the photon is injected into the unstable orbit with ra- 
dius XmiO- Conversely, photons with a large approach 
parameter are weakly deflected. Of course, the relation 
between e and i5 ensures that the SDL can be equivalently 
stated in terms of either of the two parameters. 

Let us introduce our gravitational lensing configura- 
tion. As said before, the observer is at radial coordinate 
Dql, at polar angle "do and azimuthal angle (j)o = tt. We 
will call optical axis the line connecting the lens and the 
observer. The source is assumed to be static at Boyer- 
Lindquist coordinates (i'LSj f^s, '/'s)- 

Our lens equations are provided by Eqs. I©-©, where 
we identify the final value of the azimuthal coordinate 
with the observer's one (0/ — (j)o = n), and the initial 
value with the source's one (pi = (ps- In these equations 
there are two radial integrals and two angular integrals. 
The radial integrals are solved using the SDL technique 
and expanding all coefficients to second order in a, as 
in Paper I. The results of this procedure are reported in 
Appendix A. Similarly, the angular integrals are solved 
to second order in a in Appendix B. 

Once all integrals are calculated, we have to solve Eqs. 
Q-© in terms of the source coordinates {(ps, IJ's), so that 
they are in the lens map form 



Us = fJ-siS,0 

ps = (t>s{.5-,C) 



(32) 



Note that the lens equation will be written in terms of 
the approach parameter 5 and the position angle through 
^. Through Eqs. (|30|l and (|26|l we can then go back to 
the coordinates in the observer sky (6*1, 6*2). 

In the following sections, we will calculate the critical 
curves and the caustics of the Kerr gravitational lens or- 
der by order. The procedure is indeed identical to that 
described in Paper I, save for the complication introduced 
by the additional parameter ^Lo. However, once we man- 
age to recast all equations in the best forms, the results 
remain very simple, so that a thorough discussion of the 
effects of spin and observer colatitude is possible. 



V. DERIVATION OF THE RELATIVISTIC 
CAUSTICS 

A. Zero-order caustics 

The first task is to recover the results for a 
Schwarzschild black hole, imposing the limit a ^ 0. 

Using the results of Appendices A and B to the zero- 
order, Eqs. I@J and © read respectively 
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where the new variable 
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tp = -21og5 + 21og[12(2 - VS)], 



(34) 



(35) 



allows us to put the equations in a very compact form. 
'0 actually coincides with the deflection induced by a 
Schwarzschild black hole with the same mass of our Kerr 
black hole. On the ground of this connection, we shall 
often refer to ip as "scalar deflection" in the following. 

The double signs coming from the angular integrals 
must be treated as follows: if the photon moves out of 
the source increasing its initial value of /i, then we must 
choose the upper signs, otherwise we will select the lower 
signs. These double signs are the relics of those present in 
Eqs. I@J and ^. For more details about their origin, the 
reader is referred to the Appendix B. m represents the 
number of inversions in the polar motion of the photon. 

Introducing the quantity 



■)po = t(-1)™ arcsin 



Mo 



(36) 



i~e 



we can easily solve Eqs. (1221)- (EU w.r.t. (ps and Ms to get 
the zero-order lens equation 



Ms 



T(-l)"Vl-e2sin(V^ + V;„), 



7r(l — m) + arctan 



^ tan ipo 



arctan 



^tan(V'-l- V-o) 



(37) 



(38) 



Since the azimuth is a coordinate with period 27r, we 
have eliminated the Sign[^] in front of rmr in Eq. H38|l . 
In the derivation of Eq. (EHJ from Eqs. (J23J and (J22J|, 



we have used the relations 
Ms 



^-^^l-e 



Mo 



1 - M? - C^ 



= =F tan(V' + V'o) (39) 



T(-l)"tani/;„ (40) 



and exploited the fact that the number of inversions in 
the polar motion is just the integer part of ('0 + V'o + 
7r/2)/7r. 

Let us understand the meaning of the zero-order lens 
equations. Eq. IpTZJI states that the photon performs sym- 
metric oscillations on the equatorial plane (recall that 



M = cos d) with amplitude y 1 — ^^ , which depends on 
the observer declination and the trajectory chosen by 
the photon (polar f = 0, equatorial i^ = ±1 or what- 
ever) . The number of oscillations depends on the scalar 
deflection t/i, which diverges when the approach param- 
eter (5^0. ^o is the initial condition of the oscillation, 
which depends on the observer declination. The double 
signs take into account the fact that the oscillations oc- 
cur in opposite ways depending on the starting sign of 

M- 

Eq. (|37|l is the azimuthal motion of the photon. It can 
be better understood when we choose equatorial photons 
with ^ = 1. Then it just reduces to </> = tt — '0, which 
states that the azimuthal shift is the scalar deflection 
minus tt, as expected in this simple case. Different values 
of ^ need to be analyzed by some spherical trigonometry, 
in order to understand the trigonometric functions in Eq. 

m- 

After the zero order lens equation is constructed, we 
can study the structure of critical curves and caustics. 
The Jacobian of the lens map, D, can be easily calculated 
from (jSZI) and 12HI)- We find 



= ±(- 



(-l)"^^L^^secV'oSinV' (41) 



^ = T{-ir^l-ecos{^ + ^o) (42) 



90s cos('0 -I- ■0o)secV'oSinV' 

d4>s _ £,sec'^{ip + tpo) 



(43) 
(44) 



^■0 l + ^2^an2(0.-HV'o) 

and using Eqs. (|31|l and (|35|l . we finally have 
d^i.dfjjs d^isd(j)s _ , ,„, sinV; 

d^ dtp dtp d^ ^^ ' yi^F' 

Since all transformations from (0,f) to {61,62) are 
non-singular (except for the points ^ = ±1), the solutions 
of the equation D — determine the critical curves. To 
zero order we have 



As expected, the critical curves correspond to values 
of the scalar deflection that are multiples of tt. Hav- 
ing introduced the most generic coordinate system for 
the black hole has not prevented us from recovering the 
Schwarzschild result. Through Eqs. ^5^, ^^ and l(^ 
we reconstruct the critical curves in the observer coordi- 
nates 



i?OL0i,fc(O = -^[l + efe] 



D0L62AO = ±^ yr^F [1 + Efc] 



(47) 



where 



ek = 216(2 --%/3) 



2g-fc. 



(48) 



is the separation of the critical curve from the shadow. 

We will refer to the integer number k as the critical 
curve (or caustic) order. Eqs. H47() describe a series of 
concentric rings, parametrized by ^, slightly larger than 
the shadow of the black hole and whose radius ^^ ( 1 -t- 
ek) exponentially decreases to the shadow radius with 
increasing critical curve order. 

The equations of the caustics are easily found introduc- 
ing Eq. gSJ into (^Hl-lISHIl and exploiting the fact that 
the number of inversions m coincides with k ii tp — kn. 
We have 



(-l)Vo, 0s = (l-A:)^. 



(49) 



As already known, the Schwarzschild caustics are 
point-like and lie on the optical axis. They are in front 
of the black hole (ms — fJ-o, with 4>s being an odd multi- 
ple of tt) for even values of k (retrolensing caustics) , and 
behind it (ms = —IJ'o, with 0^ being an even multiple of 
tt) for odd k (standard lensing). 

The SDL description is limited to large deflections 
(■0 > tt), thus working better and better for higher order 
caustics la, [13 • It cannot be applied to the first or- 
der one (fc = 1) whose full description can be derived in 
the weak deflection limit for sources sufficiently far from 
the lens. In what follows, we focus on caustics of order 
k > 2 and investigate how their structure is affected by 
the concomitant action of the lens spin and the observer 
declination. 



B. First-order caustics 

We now introduce first order corrections to the zero- 
order solutions found in the previous section. Starting 
from the results of Appendix A and B, we solve the lens 
equations perturbatively adding the first order terms to 
Eqs. lI2I|)-(|SHll, obtaining 



M. = Ti-ir^ji-esHtp+tp,) 



■ipk = kw. 



(46) 



T(-l)^ 



2aC 
3V3 



1 — ^2 cos 0^0 sin?/', (50) 



Second-Order Caustics 



(ps = (1 — ?7i)7r + arctan ^ tan -00 
|tan(V' + '0o) 
V' + 3%/31og(2%/3-3) 



arctan 
4a 



3^3 

1—^^ cos(V' + V'o)sinV'COsV'o 

The Jacobian of the lens equation to first order is 



(51) 



D = T{- 



_sinV^ ( 2aiy/l-^il \ 

^ ^ yr^e^l V3 ) 



(52) 



which is always solved by Eq. (|46|l . thus implying that 
the scalar deflection -0 and consequently the approach 
parameter 6 are not affected by lens spinning to the first 
order. Anyway, due to the spin dependence in Eq. (|30|l . 
first order corrections modify the separation of the criti- 
cal curves from the shadow. They read 



DoLe 



l,k 



3%/3e 



(l + efc) + 



[l + e' + e,(l-e')] V/Wl7, 



V3(i + e^)yr^ 



i?Vl^e^(l"efc)VW^ 



(53) 



where efc is still the zero-order separation defined in Eq. 

(EHl. 

Coming to the caustics, from Eqs. (|50|I - H51|I and 146() 
we get 



l^s = {-lf 


Mo, 


(54) 


0, = 7r(l - fc) - A0fc 


(55) 


A(/)fe = 4o 


[;^-log(2V3-3)] 


(56) 



So, caustics are still point-like but the alignment with 
the optical axis is now missing because of first order cor- 
rections, as already pointed out in Paper I. The azimuthal 
shift is proportional to the caustic order, it does not de- 
pend on the observer declination and is negative, thus 
implying a clockwise drift, if we look at the black hole 
from the north pole. This means that, as k is still the 
number of inversion points, prograde (retrograde) light 
rays, emitted by a source on a caustic point, perform 
more (respectively less) than (fc — l)/2 loops. Moreover, 
as the caustics drift from the optical axis and from each 
other, perfect alignment of observer, lens and source is 
not required for the enhancement of the images which 
are enhanced one at a time, as sources cannot cross more 
than one caustic point at the same time. For numerical 
values of the shift see Paper I, Table 1. 



In this section we investigate the effects of second order 
corrections in the black hole spin on the critical curves 
and caustics. Following the same steps as in the previous 
subsection, we can add the second order terms a^S^s 
to Eqs. 1)50(1 and a^(5(/)i to H51|l. Since they have quite 
long expressions, we report them in Appendix C and pro- 
ceed with the analysis of the second order lens equation. 
In fact, although the general second order lens equation 
is quite involved, it is easy to solve the Jacobian deter- 
minant equation _D = in terms of the second order 
perturbation of 0, starting from the zero order solution 
(P|l. We get 

'ipi, = kn + a Sip, 



where 



5ip 



1 
18 



(57) 



9cki3-2^ll-'si^) + 32{l-e) (58) 



and 



Cfc = -(5fc7r-t-8\/3-36). 
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(59) 



Using Eqs. (^5)1 and P?1|l we can calculate the second 
order corrections to the approach parameter S and the 
separation from the shadow e. After that, by Eqs. H26|l . 
we can derive the second order corrections to the critical 
curves given in Eq. ((53|l 



D 



OL 



h,k 



(2) 



= a2^ 

3\/3 



+efe 29- 



-2^2_2a 



^OL02^k = 



'2\ 



(60) 



+efc 



(21 - 32e2A 



1 + 2^.2 - 2e 

32^2^ 



where the zero-order separation e^ is always given by Eq. 

Plugging Eq. 157|) into the lens map, we get the caus- 
tics parametric equations up to the second order in a: 

Ai, = (-l)^^o ± c?ck{\ - Mo')'/'(l - ef'\ (61) 

0, == (1 - k)-n - A0fc - a^CfeCVl - Mo^. (62) 

As explained in Section IV Al the double sign in Eq. 1(61(1 
allows for the possibilities that the photon starts its jour- 
ney by increasing /j, or by decreasing ^, respectively, it is 
necessary to take both possibilities into account in order 
to cover the whole caustic. In agreement with Paper I 
and other works where the same results are found numer- 
ically (e.g.|l3), we get extended caustics whose shape is 
a 4-cusped astroid, with cusps in ^ = ±1 and ^ = (for 
different signs of initial /i). The extension of the caus- 
tics along /x and along (p is different. However, choosing 
appropriate coordinates centered on the caustic, it is pos- 
sible to show that the extension in the sky as seen by the 
black hole is the same along both axes (see next subsec- 
tion) . 



D. Observables related to critical curves and 
caustics 

After second order corrections to critical curves and 
caustics have been derived, we can discuss their depen- 
dence on a and do- 

First we note that the critical curves obtained adding 
Eq. (|^ to 1(351) satisfy the elhpse equation 



n,k 









(63) 



with the same origin shift as the shadow (Eq. (|22() ') and 
semiaxes given by 



^2 4-efc(4-9cfcMg) 

9 
+27cfcefc(l + efc)(3-2/x2)] 



(64) 



12\/3(l + efc) 



• (65) 



The critical curves tend to coincide with the shadow in 
the limit fc ^ 00, which corresponds to photons winding 
an infinite number of times, thus tending to the unstable 
photon orbit. The ellipticity of the critical curves is 



a^(l-M^) 



4(3 + e^.)+81cfcefc(l + £fc) 

54(1 + 6^)2 



(66) 



which is slightly higher than that of the shadow for 
the lower order critical curves, but tends to that of the 
shadow as fc ^ cx). In particular, we see that shift and 
ellipticity of the critical curves still depend on the com- 
bination asinz?o, as for the shadow. So, even the obser- 
vation of several critical curves cannot help to determine 
a and do separately. 

Let us come to the caustics. Here the situation is more 
subtle and needs to be investigated with grain of salt. 

Suppose we have no independent knowledge of the di- 
rection of the black hole spin or, at least, the direction 
of the spin is not known to any great accuracy. Then, 
the observer will construct his coordinates allowing for a 
non- vanishing position angle v for the spin axis. The un- 
certainty in V will be determinant in the following discus- 
sion. Let us thus introduce (x, -)?, cf)) as observer-oriented 
coordinates, still centered at the black hole, but with 
the polar axis perpendicular to the optical axis and the 
azimuth (p taken from the direction opposite to the ob- 
server. In general, if the observer ignores the spin axis, 
the spin axis of the black hole would have a position an- 
gle u from the polar axis as fixed by the observer. The 



n 

7 


y -^ A \ 


/ 


y ~^ "^ \ 


_v / 


1} 
0-^ 


a^^^ 



FIG. 3: Observer-oriented coordinates (1?, (^) introduced in 
the text. L is the black hole with spin a. O is the observer 
and P is a generic point, tfo is the inclination of the spin on 
the line of sight, v is the position angle of the spin. 



coordinate transformation from (/i, 0) to (i?, (/)) is 
1? = arccos [i\/\ — [i^cosv -I- [io\J'^ — /x^ cos cost' 

(67) 



-yl — /i^ sin(/)sinz^ 
arctan 



I \J\ — /i2 sin cos V 
[io\l^ — M^ cos sin z^ — /iy 1 ~ Mo^ini^) 



1 - [i^\J\ - filcOSCJ)- ^i^io 



(68) 



Fig. OiUustrates the geometrical meaning of these co- 
ordinates. 

Transforming the caustics H61|) - (|62|l from the spin- 
oriented coordinates (/i, 0) to the observer-oriented co- 
ordinates (i9, 0), and expanding to second order in a, we 
get 



ds 



^-(-l)'=A0fcyi^M!sini. 



,1 



-(-l)^^A02feMoVW^ 



-Rk 



cosz/ 



(-l)''^C-^sini/±(l-0 ' cosj/ 



(69) 



= (1 - fc)7r - A^fc^l -/i^cosi 



-2^<^feAiov/l-MoShri^ 



-Rk 



S;cosv±{-lY{l-iYi^smv 



(70) 



where 



Rk = a'ck{l-^il) = -a'{l-^il){bk^: + %V?>~iQ) (71) 
y 

is the semi-amplitude of the caustic. In fact, we can 
appreciate that, in observer-oriented coordinates, the ex- 
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tension of the caustic is the same in both polar and az- 
imuthal directions, as anticipated before for any coordi- 
nate system centered on the caustic. So, the extension is 
quadratic in the spin and is maximal for equatorial ob- 
servers, while the astroid shrinks to a single point when 
the observer lies on the spin axis. The caustic extension 
also increases linearly with the caustic order k. 

Then, we note that the angular shift of the center of 
the caustic from the optical axis is 



A, 



arccosfsini^cost 



^^l 



4av/r 



M^ 



kn 
373 



A^fcVl 



log(2V3 - 3) 



(72) 



It is linear in the black hole spin and the caustic order. 
Similarly to the semi-amplitude, also the shift is maximal 
for equatorial observers and vanishes for polar observers, 
when the axial symmetry is restored. 

The shift and the semi-amplitude of the caustics are 
very easy quantities to determine in case of observation 
of the relativistic images generated by a source cross- 
ing a relativistic caustic. In fact, if the observer is 
able to identify the source and follow its direct image 
throughout the duration of the caustic crossing event, 
then he would immediately determine the position of 
the caustic and estimate its extension. Unfortunately, 
even in these two quantities, the black hole spin and the 
observer declination always appear in the combination 
tti/l — Hg = asin^o, making the breaking of the degen- 
eracy between these two parameters impossible. On the 
other hand, it is easy to determine the order k of the 
caustic involved in the lensing event, since the ratio 



Rk 



8[fc7r-F3%/31og(2%/3-3) 
3 (SfcTT + 8V3 - 36) 



(73) 



only depends on k and increases monotonically in fc, with- 
out degeneracy between any two values. 

One possibility for the separate determination of a and 
fio arises in case the spin position angle v is known to a 
very good accuracy from independent measures. Then we 
can move to a more convenient coordinate frame where 
1/ = 0. If this is possible, looking at Eqs. (|^ and ifTUIl 
we see that the shift in the azimuthal direction is linear in 
a, while a residual quadratic shift is present in the polar 
direction, which amounts to 



1 



A0|/x<, ^1 - fil 



VoVl 



Mo 



kn 
3%/3 



log(2\/3-3) 



(74) 



Then, if one is able to measure this residual shift, one 
can extract the observer colatitude -do as 



COt^o = —J. 

At 



(75) 



Once the observer position relative to the spin axis 
is known, we can use either A^ or Rk to extract the 




FIG. 4: A typical caustic in Kerr lensing. The extension is 
the same in both directions. Having chosen coordinates such 
that the position angle of the spin vanishes, the caustic has 
an azimuthal shift A^ and a vertical shift Sk w.r.t. the line of 
sight. 



black hole spin a. However, as for the case of the direct 
determination of a from the measure of the minor semi- 
axis of the shadow, this is a higher order measure, which 
requires very accurate independent information. 

Fig. 2] shows a caustic and illustrates the meaning 
of the semi-amplitude Rk, the horizontal shift A^ and 
the vertical shift Sk- The picture is done for a standard 
lensing caustic (fc odd) with -do > n/2, so that the caustic 
is displaced upward (see Eq. (|H^ ). 

As usual we can trust our results as long as the per- 
turbative terms remain small. In extremal or close-to- 
extremal Kerr black holes, higher orders in a would play 
a major role in the critical curves and caustics profile. In 
that case, the degeneracy between a and z9o can be prob- 
ably broken also through the determination of the exten- 
sion and position of the caustics or through the analysis 
of the critical curves. However, in the literature there is 
no investigation of Kerr black holes with high spin that is 
deep enough to allow a comparison with our perturbative 
results for low spins. 



VI. GRAVITATIONAL LENSING NEAR 
CAUSTICS 

A. Position of the relativistic images 

Although in our picture the images cannot be found 
analytically for arbitrary source positions using the lens 
mapping that we have derived, they can be actually found 
for sources in the neighbourhood of a caustic. This is in- 
deed the most interesting case, as the relativistic images 
are highly magnified and become observable only if this 
event occurs. Assuming that the angular distance be- 
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tween the source and a caustic of order k is of the order 
of 0^ (thus comparable with the caustic semiaxis), we can 
write the source position as 



Ms 



(-l)Vo + a^<5Ms, 



= (1 - k)TT - A(j)k + a^ 



(76) 



(77) 



In this assumption, the images will be very close to the 
critical curve of order k. Then the scalar deflection will 
be 



ip = kn + a Si/j. 



(78) 



Plugging the last equation into the lens map written 
up to corrections of second order in a and inverting with 
respect to S^s and Scfig, we get 

1 



+{5kTr + 8\/3 - 20)(1 - (1 - fio^)^^)] , (79) 



e 



92-24\/3- l5kTT-9S^p 



9^1 -Mo" 

+2^o"(5/c7r + 8\/3-36) 

+ (5fc7r + 8\/3 - 20)(1 - Mo")C^ 



(80) 



Solving (|Hn|l with respect to Sip and plugging its ex- 
pression into (|79|l . we find 



6fis^ = -S{-lf{l - /io')(<50s + CfeVT^M^Ov/^^I^, 

(81) 
where Cfc is given by Eq. (|5^ and S = =f(— I)''- This 
equation can be more conveniently written in terms of 
observer-oriented coordinates {"& sift's)- Supposing that 
the position angle of the spin has been well established by 
observations of the shadow or by the shift of the caustic 
itself, we put ly ^ for simplicity and write 



where Rk is the semi-amplitude of the caustic given by 
Eq. H71|l . The solutions of this equation for arbitrary 
source positions ((5i?s, (50^) determine the relativistic im- 
ages generated by the Kerr black hole. As the roots of 
Eq. l|S^ are found squaring both its sides, the solutions 
of the squared equation satisfy the original one only for 
one choice of S. S is directly related to the half-sky 
where the image appears. In fact, we recall that the pa- 
rameterization (|26|l has an ambiguity in the sign of 62. 
This ambiguity can be solved observing that the photon 
reaches the observer from the upper side of the black hole 
if S is positive and from the lower side if S is negative. 
This fact can be easily established remembering that in 
all our equations the upper signs hold when the photon 
leaves the source by increasing its /i coordinate. Then, 
if its polar motion undergoes one inversion (k — 1), the 
photon reaches the observer from above and we coher- 
ently have S = 1. On the other hand, if the lower signs 
hold, the photon begins its motion decreasing its fi coor- 
dinate. With one inversion, it reaches the observer from 
below and coherently we have 5" = — 1. The same rea- 
soning can be repeated with an arbitrary number k of 
inversions in the polar motion. 

It can be easily verified that Eq. (|86|l has four real so- 
lutions if the source is inside the caustic and only two real 
solutions if the source is outside. Once the coordinate £, 
(which, we recall, represents the cosine of the position 
angle) of the image is known, Eq. H80|l can be used to 
determines the value of Sip (perturbation of the scalar de- 
flection). However, it is important to stress that Eq. H86() 
determines ^ to zero order only. Therefore, though the 
positions of the images in the observer sky are generically 
given by 



DoL01 



3^3 



e(l + efc) 



£|{l0-4e"(l~M^)-4M^ + 6,[27<5^ 
+58- 64^(1 -fil)-l6fil]}, (87) 



^s^^ -i-iySk + 6ds 



4>s ^ (l - k)TT - A4>k + S4>s, 



with 



S'ds = —a 



2 ^Ms 



Vi^m! 



S4>s = a^S(l)sV^ - Mo- 



(82) 
(83) 

(84) 
(85) 



Then, we can write Eq. (|81|l directly in terms of these 
coordinates as 



S^s^^Si-l)HS$s + RkOV^^^, (86) 



,3\/3 



DoL02= S^Vl-e{l + ek) 



-a5eVl-e'v/T^M!(l - ck) 

-s£=VT^H2-4e{l-^^l) + 4^^l 

+ek [27S^ + 42 - 64(1 - /lo)^'] } , 

to the second order in a, only a zero order expression of 
^ is actually available. So, the position of the images is 
accurate only to zero order in a and is given by 



DoL0i^ -^m + ek) 
,3\/3 



(89) 



Dol92^ S^Vl^{l + ek). (90) 
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To zero order, we see that the images of order k he 
along the critical curve of order k (we remind that e^ is 
just the separation of the critical curve of order k from 
the shadow H48|l ). with position angle determined by the 
solutions of Eq. H86|l . If a more accurate theoretical 
prediction of the images position (including first order 
corrections) is needed, it is necessary to push the lens 
equation to the third order. Indeed this would be a wor- 
thy (though heavy) task since the equation for the images 
(|5S|l depends on a only through it!^. As noticed before, 
this quantity only depends on the projection of the spin 
on the line of sight. So, once more, the observables (in 
this case the positions of the images) only depend on 
asindo to the lowest order. However, contrarily to the 
former observables, the positions of the images could be 
detected to an accuracy sufficiently high to be sensitive 
at least to first order corrections in a. So, it would be 
indeed desirable to check whether the positions of the 
images may help to break the degeneracy between the 
absolute value of the spin and its inclination on the op- 
tical axis. 



B. Magnification 

The magnification is defined as the ratio of the an- 
gular area of the image and the corresponding angu- 
lar area of the source. The angular area of the image 
is simply \d9idd2\, while the angular area of the source 
is \sin'dsd(j)sd'ds\ or Isin-d sd(l)sd'ds\ if one uses observer- 
oriented coordinates. Then the magnification can be cal- 
culated as |sin^s|~^ times the inverse of the Jacobian 
determinant of the lens application in the form 






(91) 



Following the same approach of Paper I, we can find 
the expression of the magnification for sources in the 
neighbourhood of caustics exploiting the available rela- 
tions (JHSJ-IHII) and IZnj-llHni) to get 



and |57|l -P ll 












(92) 



(93) 



Then the perturbation of the scalar defiection 5%l) and the 
cosine of the position angle ^ play the role of intermediate 
variables between the source coordinates {'dg , ips ) and the 
image coordinates in the observer sky [9i, 62). 

Since smd sd<t)s — d{64>s) and d^s — d{S-ds) to the low- 
est order, the Jacobian of the map (|91|l reduces to 



d{9i,02) d{S^,0 



d{9i,92) 



d{S4>,0 



(94) 



where we have used the matrix notation 



diy,,y2) _ ^ ^ 



9(xi,a;2) 



Oyi dy2 
dxi dx2 



(95) 



As the derivatives and the Jacobian matrix have very 
involved expressions, we do not go too much into detail 
and only report here the two eigenvalues of the Jacobian 
matrix 



Ar 



2Doi 
3V3( 



efc 



2DolDo 



where 



27V3(1 

Do= (-l)'^y{9cfe 
+32(1 - e^) 4 



■Cfc 



i-2^il- 3|2 



(96) 



(97) 



18(5iA| 



In a first approximation A^ only depends on the caustic 
order k and is always positive. On the other hand At 
vanishes at caustic crossing (see. Ea. (|5H|l l Following 
Paper I, we will call A^ and At, respectively, radial and 
tangential eigenvalues, although they are such only in the 
limit a — > 0. Taking into account that the flux received 
by the observer is Dj^g/DQg times the flux received by 
the black hole, the radial and tangential magnifications 
are 



D 



^ir 



fJ-t 



OS 



DlS Ar 



Dos 1 

Dls I At I 



(99) 



(100) 



while the total magnification is given hy fj, — ^r-Mt- 

An interesting thing to note is that the radial magni- 
fication is completely independent of a and ^o- It is just 
the same as in the Schwarzchild black hole case. On the 
other hand, the tangential magnification is sensitive to 
the caustic structure, which can be seen more clearly if 
we plug the solution of the lens equation H80() for 6ijj into 
Eq. (jnU- In fact, we have 



M* - (-1) 



k^OS 



3^3(1 + efc)e 



Dls 2Dol (Rk^' 



(101) 



where the (—I)'' accounts for the parity of the image 
and ^ must be determined solving Eq. (|86|l . The whole 
dependence of the magnification on the black hole spin 
and the observer declination is through the caustic semi- 
amplitude Rk, where they appear in the usual combina- 
tion asmiSr,. 
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C. Relativistic images around Sgr A* 

In this subsection we want to complement the discus- 
sion about the detectabihty of relativistic images done in 
Paper I by some additional considerations. Indeed there 
are many factors that contrast the positive detection of 
relativistic images around Sgr A*. The photons with the 
right incident direction for performing a complete loop 
around a black hole and then reach the observer are very 
few, because a slight perturbation in the incident trajec- 
tory results in a very different outgoing direction. More- 
over, during their journey, photons may be scattered or 
absorbed by the accreting matter surrounding the super- 
massive black hole. Finally, the photons surviving up to 
the observer must be recognized and distinguished from 
the noise coming from the environment. 

Scattering and absorption from accreting matter are 
strongly model-dependent and cannot be easily quan- 
tified without non-trivial assumptions on the infalling 
plasma physics. We are not going to face this problem 
here, since it demands an extensive investigation beyond 
the purpose of this work. 

On the other hand, our gravitational lensing analy- 
sis allows us to give sharp answers on the brightness 
and spatial properties of the images. In Paper I, we 
have suggested that the observed Low-Mass X-ray Bi- 
naries (LMXB) orbiting around Sgr A* provide an ideal 
population of sources for the gravitational lensing in the 
SDL ■23. Of course we need to resolve the shadow of 
Sgr A* in order to identify relativistic images around it. 
This requires a resolution of the order of the /ias, which 
is just one step beyond the limit reached in the radio 
band. In the X-ray band, projects of space interferom- 
etry which could reach res olutions even better th an ^as 
are under study (MAXIM, 'http://maxim.gsfc.nasa.gov). 
When such projects will become reality, a complete imag- 
ing of Sgr A* will be possible and the relativistic images 
could be distinguished. 

Apart from spatial resolution, which can be attained by 
realistic future projects, in order to detect a signal in the 
X-ray band from a relativistic image, we need a sufficient 
energy flux. With an intrinsic luminosity Ls ~ 2 x 10'^'^ 
ergs s~^ in the band 2 — 10 keV, emitted by a surface with 
radius Rs = 100km, LMXBs are as powerful sources as 
Sgr A* itself but enjoy a much higher surface brightness 
[23 . If one of these sources crosses a relativistic caustic 
of order fc, the angular area of the resulting relativistic 
image is the original source area irR^/DQg multiplied 
by the magnification factor /x. As long as the source is 
inside the caustic, the magnification stays higher than a 
minimum value corresponding to a source located at the 
center of the caustic. The central magnification has been 
calculated in Paper I and amounts to 



D2 



27efc(l + efc) 



d2 n2 



^Rk 



(102) 



source is inside the caustic. 

For a detector with collecting area Ao, the observed 
flux, taking into account an absorption factor e = 0.158, 
deduced from Ref. [2g, is thus 



Fk=e 



Ls 

AttRI 






Ai 



(103) 



With DoL = 8kpc, Mbh = 4.3 x 10*^ Mq |2J and 
Dos — DoL (since Dls < Dql), we have 

...2..io-.ws-.(ji^)"'(5;i,)-X^) 

(104) 
for a source crossing the caustic of order fc = 2 and a black 
hole spin a — 0.02 [23. This flux is independent of the 
source radius, as long as the source is much smaller than 
the caustic extension, as in our case. We have considered 
a collecting area Ad — lOOm^ which might be realisti- 
cally obtained by future space detectors. The count rate 
for photons in the considered band (with average energy 
6 keV) is thus of the order of 2.4 x 10~^ s~^, which is 
comparable to the counts usually reported as positive de- 
tections by the Chandra satellite for faint sources |23 . l2a |. 
Of course, such a high value for the count rate can only 
be achieved with a collecting area as large as that we 
have considered here, which is roughly 100 times larger 
than that of Chandra. 

Sgr A* itself emits in the X-rays and provides a back- 
ground noise to the signal of a relativistic image. The 
image of an LMXB is entirely contained within a single 
pixel of a hypothetic detector where every pixel covers 
Ifiasxlfias of sky. We can estimate the noise due to Sgr 
A* considering that its intrinsic luminosity is of the same 
order as Ls 23, 26], but its emitting region has a radius 
Rsgr of the order of 100 Schwarzschild radii. Then, every 
pixel is affected by a noise from Sgr A* of the order of 



Fsgr — £ 



Ls 2 



uj^Ad, 



where Up is the size of the pixel. We thus have 
^.,..3.7xl0-ergss-(^)^(^) 



(105) 



(106) 



for each of the four relativistic images present when the 



which is roughly 600 times smaller than ^2. This proves 
that the background from Sgr A* is indeed negligible for 
relativistic images of order 2 if one has sufficient resolving 
power. It is also important to stress that these estimates 
has been calculated considering the minimum magnifica- 
tion fXc for a source inside a caustic. When the source is 
close to a fold or a cusp, the brightness of the relativistic 
image can be sensibly higher. 

We conclude this discussion mentioning that the 
brightness of relativistic images of order 3 is 0.016i^2, 
which allows a marginal detection w.r.t. the noise by Sgr 
A*, while relativistic images of higher order are too faint 
to be detected, at least for the configuration examined 
here. 
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VII. CONCLUSIONS 

This paper completes the cycle of papers devoted to 
the study of gravitational lensing by Kerr black holes in 
the Strong Deflection Limit. After the first pioneering 
work of Ref. [l5l | , where equatorial lensing was reduced 
to the same problem already solved for spherically sym- 
metric black holes j5|, in Ref. [T^ we managed to make 
a complete analytical treatment of Kerr lensing for equa- 
torial observers, introducing a perturbative expansion in 
the spin a. In this work we have extended that idea to 
Kerr lensing with a generic observer. Though the strat- 
egy is essentially unchanged, the introduction of a new 
parameter (the inclination of the spin or equivalently the 
observer colatitude do) has increased the difficulty of the 
derivation. Nevertheless, our investigation has reached 
its objective: a basically simple and accurate description 
of Kerr lensing phenomenology with arbitrary observer 
position. 

An essential summary of the main results obtained in- 
cludes: the shape of the shadow of the black hole (|21|l : 
the shape of all critical curves (|63|l : the shape and po- 
sition of the caustics (Eqs. (|55|l and (|7n|) '): the position 
of the images (Eqs. (jSHJ-lEnjl with Eq. ^) and their 
magnification (Eqs. H99|l with (|96|l and HlUl|l 'l for sources 
close to a caustic. 

To the second order in a, the shadow of the black hole 
and the critical curves are ellipses slightly displaced from 
the black hole position. The ellipticity is slightly higher 
in critical curves than in the shadow. The caustics are 
displaced from the optical axis and show the character- 
istic 4-cusped astroid shape with the same extension in 
both directions. The caustic shrinks back to a single 
point when the observer lies on the spin axis, restoring 
the axial symmetry. There are two additional images 
when the source is inside a caustic. 

The fundamental fact that emerges is that all observ- 
ables to the lowest order are functions of a sin 'do , which 
represents the projection of the black hole spin on a plane 
orthogonal to the line of sight. These observables include: 
the shift and the ellipticities of the shadow and of crit- 
ical curves; the shift and the extension of the caustics; 
the position and the magnification of the images. 

The degeneracy between the absolute value of the spin 
and its inclination on the line of sight can only be broken 
by next-to-leading order terms in all observables. This 
has been explicitly shown considering the shadow and 
critical curves semi- axes and the caustic vertical shift. 
These are second order contributions to zero-order quan- 
tities, thus requiring extremely accurate measures, which 
may be very challenging. For example, if the black hole 
spin is a = 0.1, in order to break the degeneracy we need 
a relative accuracy of order a^ = 0.01 in the measures. 

The most promising way to break the degeneracy is 
through higher order corrections to the positions of the 
images. In fact, our second order treatment is only suffi- 
cient to determine the position angle of the images to zero 
order in a. Indeed the first order corrections are likely 



to be at reach of future VLBI observations, but unfor- 
tunately they require at least a third order treatment of 
Kerr lensing in order to be determined analytically. This 
could represent the main target for future theoretical de- 
velopments of our methodology. 

Of course, if the black hole spin is close to the extremal 
value a — 0.5, the degeneracy breaking terms arising from 
higher orders in a grow to the same size as the lowest 
order contributions and the problem would not be the 
degeneracy between a and do but the correct theoretical 
interpretation of the observations in a non-perturbative 
frame, in order to perform a safe parameters extraction. 
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APPENDIX A: RESOLUTION OF RADIAL 
INTEGRALS 

This appendix reports the calculation of the radial in- 
tegrals appearing in the geodesies equations Q) and |(SJ). 
The double signs remind us that the integration along 
the whole trajectory of the photon must be performed 
in such a way that all pieces bounded by two consecu- 
tive inversion points must sum up with the same sign [23| . 
Gravitational lensing trajectories have only one inversion 
point in xq, the closest approach distance. Thus we just 
have to sum the contributions due to two branches (ap- 
proach and departure) . These two branches of the photon 
trajectory are actually related by the time-reversal sym- 
metry, so that the results of the whole radial integrals 
are just twice the contributions covering the departure 
branch. Summing up, the radial integrals reduce to 



J y/R 

Xq 

oo 

J AVfl 



(Al) 
(A2) 



where we have neglected the corrections due to the finite- 
ness of DoL and D^s, thus extending the integration do- 
main to -|-(xi. The resolution by the SDL technique can 
be read from the appendix A of Paper I, since the only 
change comes when we replace J and Q by their new ex- 
pressions containing fio- Thus we can directly jump to 
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the results, which read 



/i = -fli log (5 + bi 
h = -a2logS + b2. 



(A3) 

(A4) 



The coefficients expanded to second order in a are 



fli 



16 



3V3 27 81V3 



a^T + Afil + be) 
8 



(A5) 



6i = ai log[12(2 - V3)] - —{dVs - 6)^2(1 - f) (^g) 



"^^71 + 3«^+9V3 



02(3 + 2^2^5^2 



(A7) 



62 =a2log[4V3(2V3-3)i+^ 



--aC"[9 - 2\/3 + 3V31og(2V3 - 3)] 
9 



+—a^ {26\/3 - 16 - 2V3 log(3) + 8fil - 12^3/^2 

+ 12(3 - fil) log(2\/3 - 3) - S,^[38V3 - 20 

+5(1 + 2\/3) log(3) + 30 log(2 - Vs)] } (AS) 



with e = evi^T^- 



APPENDIX B: RESOLUTION OF ANGULAR 
INTEGRALS 

This appendix is devoted to the resolution of the an- 
gular integrals 



J2=± 



CSC^l? 



d^. 



(Bl) 

(B2) 



Introducing the variable /i — cos??, the two integrals 
become 



Ji=± 

J2=± 
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/a 



id/^ 



zd^, 



where 



(i-M2)v/e; 
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(B3) 
(B4) 

(B5) 

(B6) 

(B7) 



and we have already replaced J and Q with J,„ and Q™, 
coherently with the fact that we only retain terms that 
are logarithmically diverging or constant in the approach 



parameter 5 (or equivalently in the separation from the 
shadow e). 

Q^ has two zeros in // = ±/i+. Then the photon per- 
forms symmetric oscillations of amplitude ^+ w.r.t. the 
equatorial plane. It is useful to write the explicit ex- 
pressions of ^+ and /i_ in terms of the spin a and the 
position parameter ^. Using Eqs. lfT7|) - ((TH|) in Eq. HB6|) 
and expanding to the second order in a, we find 



/i+ = 



1-^2 

3\/3 



l + aA 
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where 



A.= 



2e(i-Mg)(i-a 

3\/3(l - ^2) 



(B8) 



(B9) 



(BIO) 



In a first approximation, the oscillation amplitude /i4 



is y 1 — ^2 , plus corrections due to the black hole spin. 
Note that the minimal amplitude of the oscillations is 
obtained for ^ = ±1, which gives ^+ — |/io|. Purely 
equatorial trajectories with /i+ = are involved in gravi- 
tational lensing only if the observer itself lies on the equa- 
torial plane. On the other hand, polar photons (^ = 0) 
perform oscillations with maximal amplitude /i+ = 1, 
touching the poles of the black hole. 

Now it is convenient to introduce a new integration 
variable z = A*//^+j which allows to eliminate the de- 
pendence on a in the integration extrema. The integrals 
become 
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zdz. 



with 
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In order to perform the angular integrals, it is wise 
to expand the integrands to second order in a and then 
integrate. The primitive functions read 
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Similarly to radial integrals, the angular integrals ap- 
pear with double signs reminding that they must be per- 
formed piece by piece between any two consecutive inver- 
sion points and all contributions must be summed with 
the same sign ^20] . The integration covers the whole tra- 
jectory of the photon, which may perform several oscilla- 
tions around the equatorial plane. The integration must 
start from the source position Zg = fig/fi-^- and must end 
at the observer position Zq = fio/fJ'+- Let us indicate by 
m the number of inversion points in the polar motion 
touched by the photon. Still we must consider two pos- 
sibilities depending on the direction taken by the photon 
starting from Zg. In fact, we may have a trajectory in 
which z is initially either growing or decreasing. In the 
first case, the first pieces of the angular integrals cover the 
domain [zg, 1]. After that, we have m — 1 integrals cov- 
ering the whole domain [—1, 1]. All these integrals must 
be taken with the same sign so that they always sum up. 
Finally, if m is even, the photon reaches Zq with growing 
z and the last piece covers the domain [— l,Zo], other- 
wise z is finally decreasing and the domain is [zo, 1]. The 
total angular integrals are thus given by the sum of all 
these contributions covering the domains just described. 
Exploiting the primitive functions (jB14|) and ljB15|l . we 
can express each integral as (in the following, i takes the 
values 1 or 2) 

J, = Fj^ (!) - Fj^ (zg) + (m - 1) [Fj^ (1) - Fj^ (-1)] 
+Fj^{z,) - Fj^{-1) (B16) 

for 771 even and 

J, - Fj^ (1) - Fj^ (zg) + (m - 1) [Fj^ (1) - Fj, (-1)] 
+Fj^{l) - Fj^izo) (B17) 

for m odd. 

Noting that both primitives are odd functions of z, we 
have Fj. (— 1) = —Fj.{l) and we can express the angular 
integrals in the compact form 

J,; = T [Fj, (zs) - i-iy^Fj^ (zo)] + 2mFj, (1). (B18) 

The (—1)™ ensures that the sign of the Zo-term is the 
same as the z,-term if the number of inversions is odd 



and is opposite if m is even. We have also introduced a 
double sign to take into account the possibility that z is 
initially decreasing from the starting value Zs- 

For future reference, we also write the explicit values 
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APPENDIX C: SECOND ORDER 
CONTRIBUTIONS TO THE LENS EQUATION 



In this appendix we report the expressions for Sfj, 



(2) 



and I 



a(2) 



which must be added to Eqs. H50() and (|51|l to 



obtain the second order lens equation. They read 



Sfii'^ =T(-ir 
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U{i^ -l)cosipi 



+ (l + 3e^)cosV'o 
+ (l-|2)cos(2V'i -l-?Ao) sin^AJ 
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arctan(^tan V'o) ~ arctan(^tanV'i) 
2 + (l-^"2)cos2?Ai 



9(cos2 -01 + ^2 sin2 t/,^) 
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864^(cos2 ipi + ^2 sin^ V'l)^ ^ 



■EP^i' 
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where ipi — tp + ijjo, ^n = 5-0 + 8\/3 — 20 and 

Po = 64 cos"^ i/)i cos'^ Vo sin -0 (C3) 

pi = 384 cos V'l + 2 (7 + 4 cos 2ipo + cos Aipo) sin 2f/'i 
-(ll + 20cos2'(/'o + 5cos4V'o)sin4V'i 
+4 [384 + (14 sin 2-00 + 5 sin 4^0) cos^ Vi] cos 20i 
+96(12 -cos^tAiCOS^Vo sin ?Ao) (C4) 

P2 = 768(1 -cos 4^1) -8(9 + cos2?Ao)cos^0oSin2V'i 
+ ( 1 3 + 28 cos 2-00 + 7 cos A^o ) sin 40i 
+8(9 — 11 cos 2^1) cos^ Vi sin 2?/;o 
+4(9 + 7cos2V'i)sin^0isin4V'o (C5) 

P3 = 384(3 + cos 4'0i) - (5 + 3 cos 4-0o) sin 4-01 
+2(9- coa^ipo + 24cos20oSin^ Vi) 
-20(sin^ 2V'i sin 20^0 + sin^ 0i sin 40'o) 
-12(128 + sin^ -01 sin40o) cos 2-0i. (C6) 



[1] C. Darwin, Proc. of the Royal Soc. of London 249, 180 [11] 
(1959). 

[2] J. P. Luminet, A&A 75, 228 (1979). 

[3] H.C. Ohanian, Amer. Jour. Phys. 55, 428 (1987). 

[4] V. Bozza, S. Capozziello, G. lovane, and G. Scarpetta, [12] 
Gen. Relativ. Gravit. 33, 1535 (2001). 

[5] V. Bozza, Phys. Rev. D 66, 103001 (2002). [13 

[6] E.F. Eiroa, G.E. Romero, and D.F. Torres, Phys. Rev. [14 
D 66, 024010 (2002); A. Bhadra, Phys. Rev. D 67, [15 

103009 (2003); E. F. Eiroa, Phys. Rev. D 71, 083010 [16 
(2005); R. Whisker, Phys. Rev. D 71, 064004 (2005); 
A. S. Majumdar and N. Mukherjee, Int. J. Mod. Phys. [17 
D 14 , 1095 (2005); J. M. Tejeiro and E. A. Larranaga, [18 

lir^/ 0505054 E.F. Eiroa, gr-qc/0511065 K. K. Nandi, [19 
Y.-Z. Zhang, and A. V. Zakharov, gr-qc/0602062j K. [20 
Sarkar and A. Bhadra, gr-qc/0602087 

[7] S. Frittelh and E.T. Newman, Phys. Rev. D 59, 124001 [21 

(1999); S. Frittelh, T.P. Khng, and E.T. Newman, Phys. 
Rev. D 61, 064021 (2000); M.P. Dabrowski and F.E. [22 

Schunck, Astrophys. J. 535, 316 (2000); K.S. Virbhadra 
and G.F.R. Ellis, Phys. Rev. D 65, 103004 (2002); V. [23 
Perlick, Phys. Rev. D 69, 064017 (2004); V. Perlick, [24 
Living Rev. Rel. 7, 9 (2004); P. Amore and S.A. Diaz, 
gr-qc/0602106 

[8] K.S. Virbhadra and G.F.R. Ellis, Phys. Rev. D 62, [25 

084003 (2000). [26 

[9] B. Carter, Phys. Rev. 174, 1559 (1968). 
[10] H. Falcke, F. Melia, and E. Agol, ApJ 528, L13 (1999). 



S.U. Viergutz, A&A 272, 355 (1993); A. de Vries, Class. 

Quant. Grav. 17, 123 (2000); R. Takahashi, Astrophys. 

J. 611, 996 (2004); K. Beckwith and C. Done, MNRAS 

359, 1217 (2005). 

C.T. Cunningham and J.M. Bardeen, ApJ 183, 237 

(1973). 

K.P. Ranch and R.D. Blandford, ApJ, 421, 46 (1994). 

W. Hasse and V. Perlick, gr-qc/0511135 

V. Bozza, Phys. Rev. D 67, 103006 (2003). 

V. Bozza, F. De Luca, G. Scarpetta, and M.Sereno, Phys. 

Rev. D 72, 083003 (2005). 

M. Sereno, MNRAS 344, 492 (2003). 

V. Bozza and L. Mancini, ApJ 627, 790 (2005). 

V. Bozza and M. Sereno, gr-qc/0603049 

S. Chandrasekhar, Mathematical Theory of Black Holes, 

Clarendon Press, Oxford (1983). 

R.H. Boyer and R.W. Lindquist, Jour, of Math. Phys. 8, 

265 (1967). 

A.F. Zakharov, A. A. Nucita, F. De Paolis, and G. In- 

grosso, astro- ph/0411511, 

M.P. Muno et al., ApJ 622, L113 (2005). 

A.M. Beloborodov et al., 'astro-ph/0601273 F. Eisen- 

hauer, R. Genzel, T. Alexander et al., ApJ 628, 

246(2005). 

S. Liu and F. Melia, ApJ 573, L23 (2002). 

F.K. Baganoff et al, ApJ 591, 891 (2003). 



